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Abstract 

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner 
(OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. 

Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) 
laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having 
steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Flz, and phase 
angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input 
velocity (L/ air , at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to 
prescribed oscillatory “peak-to-peak” velocity inputs superimposed on steady inputs. 

Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Flz Hot- 
Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 
Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, 
and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. 
The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 
1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 
1 00 d(L/ air / L/ a j r 0 Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude 
(and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. 

For the C 2 H 4 /N 2 -air system, from 8 to 20 Hz, DFW is constant at 8.52 + 0.20 (% weakening)/Pa. This 
reflects a “quasi-steady flame response” to an acoustically induced dl/ air /dP. Also, it is surprisingly 
independent of C 2 H 4 /N 2 mole fraction due to normalization by SS FS. From 20 to ~ 150 Hz, the 
C 2 H 4 /N 2 -air flames weakened progressively less, with an inflection at ~ 70 Hz, and became 
asymptotically insensitive (DFW ~ 0) at ~ 300 Hz, which continued to 1600 Hz. The DFW of CH 4 -air 
flames followed a similar pattern, but showed much greater weakening than C 2 H 4 /N 2 -air flames; i.e., 
the quasi-steady DFW (8 to ~ 15 Hz) was 44.3 %/Pa, or - 5x larger, even though the 0 Hz (SS) FS was 
only 3.0 x smaller. The quasi-steady DFW's of C 3 H 8 -air and C 2 H 6 -air were intermediate at 34.8 and 
20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but 
non-linearly with respective SS FS's. Notably, the DFW profile for C 3 H 8 -air fell more rapidly in the 
range > 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical 
kinetics, and/or 0 2 transport to a flame that moved closer to the fuel-side. 
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In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and 
correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For 
reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 
Flz) than they did at low frequencies (< 15 Hz), where quasi-steady weakening appears to dominate. 
Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just 
exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed 
mechanistic understanding is needed in the fall-off region. 


Introduction 

Turbulent non-premixed hydrocarbon-air flames are frequently exposed to strong acoustic fields in 
practical combustion devices, but research in characterizing the fundamentals of flame perturbation has 
been sporadic [1-5]. Although our applied objective is to understand key physical/chemical 
flameholding processes in hydrocarbon-fueled airbreathing scramjets (see [22] for a detailed review of 
flameholding processes and ground-test effects), such engines represent but one of many different 
applications that may be affected. Despite the surge of interest in unsteady flames during the early 
1990's [3-15], much remains to be understood regarding specific physical/chemical effects of acoustic 
oscillations on the structure and extinction of even the simplest dynamically strained diffusion flames. 
Most of the current knowledge stems from large-activation-energy asymptotic analyses [6,8,9] and 
numerical simulations [5,7,11,12]. The few known experimental studies have emphasized detailed non- 
intrusive measurements of species, temperature, and velocity, and resultant axial strain rates in 
methane-air and propane-air systems at relatively low frequencies [10,13-15]. 

This paper presents an idealized experimental study of four gaseous hydrocarbon-air systems, using 
an Oscillatory-input convergent-nozzle Opposed Jet Burner (OOJB) system developed from earlier 
steady-state studies [1 6-21 ]. It features a substantially improved analysis of the ethylene-air and 
methane-air dynamic extinction data presented in our previous paper, AIAA-2003-4634 [26], and it also 
contains new sets of ethane-air and propane-air data that exhibit intermediate behavior. Both features 
lead to a simpler and broader set of conclusions, compared to our initial publication. 

Thus, about 600 dynamic-extinction “Flame Strength” measurements were obtained on unanchored 
(free-floating) laminar Counterflow Diffusion Flames (CFDF's), stabilized by steady input flows, and 
perturbed by superimposed sinusoidal velocity inputs (mostly in-phase). We define steady-state Flame 
Strength, FS, as the maximum average air input velocity (L/ air , at nozzle exit) a CFDF can sustain before 
it extinguishes suddenly, due to increasing net heat loss and a precipitous decline in temperature [21]. 
Although the steady-state strength of a flame can, in principle, be defined more closely by measuring 
the maximum steady axial strain rate near the airside edge of the flame [21 ], such a measurement is 
inadequate for dynamic Flame Strength. This is because, under dynamic conditions, the input velocity 
gradient (strain rate) field varies complexly with axial position, and with the magnitude and frequency of 
applied input velocity oscillations (pk/pk L/ air ). 

Although Dynamic FS might be gauged by a sinusoidal-maximum axial strain rate near the airside edge 
of a CFDF, that would only indirectly reflect maximum radial strain rate in the all-important flame core. 
And because the entire transient strain field varies complexly with the magnitude of sinusoidal input 
(pk/pk U a \ r ) and f, definitive flow-field measurements and deductions of transient temperature and trans- 
port effects are difficult, and unavailable. Finally, although pk/pk L/ air input magnitudes are theoretically 
derivable from measured acoustic pressure magnitudes (pk/pk P) when local acoustic impedance is 
accounted-for, such complex flow-field characterizations have not been done. Instead, global extinction 
parameters and Probe Microphone pk/pk P cold-flow measurements were used to assess Dynamic FS. 

Thus measurements of dynamic extinction were supported initially by limited HotWire, and later by 
extensive Probe Microphone, cold-input-flow calibrations. The latter enabled characterizations of the 
dynamic weakening of each hydrocarbon-air system at frequencies from 4 to 1 600 Hz. Quasi-steady 
weakening is measured at low frequencies (e.g., <_15 Hz), where the maximum L/ air attained slightly 
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exceeds the l/ air for steady-state extinction. At progressively higher frequencies, a transition (decrease) 
in dynamic flame weakening occurs, which is known to be accompanied by increasing phase lags in 
diffusive transport. Finally, at still higher frequencies, flames regain their original “steady-state strength” 
as they become unresponsive to the rate-limited molecular diffusion of flame reactants and products. 

The present data analyses and comparisons with published analytic and numerical results offer new 
insight on: The dynamic strain rate effects on quasi-steady extinction; the apparent effects of diffusive- 
transport phase lags with increasing frequency, in which diffusion flames might weaken even further as 
they depart from the quasi-steady state, and then “re-strengthen” at higher frequencies; and, the 
eventual development of complete insensitivity at high-frequencies. 

Finally, with particular application to airbreathing scramjet engines, we recognize that critical early 
stages of localized subsonic “incipient” flameholding are typically associated with fuel injection 
processes that “feed” reacting recirculation zones having “sufficient” residence-time distributions [22]. 
Such incipient flames may be weakened and extinguished by acoustically-perturbed, frequency- 
sensitive diffusion processes, such as Fl-atom diffusion in a “laminar-flamelet-like” reaction zone. 
Because such limitations may cause the loss of “robust” flameholding, and possibly generalized 
flameout, it is important to characterize these dynamic effects and begin to assess the potential for their 
occurrence in scramjet designs. Thus it appears that detailed measurements of low frequency (< 200 
Hz) acoustic fields are also needed to assess possible effects on critical flameholding in ground-based 
and flight tests of scramjet engine configurations. 


Experimental Approach 

The experimental approach follows earlier extensive studies of steady-state CFDF's using various 
nozzle and tube-OJB's [21], an initial study of the acoustic weakening of CFDF's using an oscillatory- 
input 7.2 mm nozzle-OJB system [26] unchanged from the present OOJB system, and a more recent 
Work-in-Progress Poster presentation featuring a significantly improved analysis of dynamic flame 
weakening [27]. The current results include four major sets of data on the dynamic extinction of 
gaseous hydrocarbon-air systems, one of which was nitrogen diluted. The reanalyzed data, obtained 
as functions of applied pk/pk voltage magnitudes to twin speaker-drivers, are respectively combined 
with Plot Wire and Probe Microphone (PM) dynamic flow calibration data obtained under essentially 
identical conditions, except the PM speaker-calibration flows were significant (5 SLPM), and of course 
cold. Thus dynamic flame weakening is quantified in terms of axially applied sinusoidal velocity or 
“equivalent” pressure oscillations at the nozzle exits. The combined results thus provide idealized 
characterizations of flame weakening effects due to forced oscillations of fuel and air input velocities on 
the subject laminar counterflow diffusion flame (CFDF) systems. 

The dynamic OOJB system shown in Fig. 1 consists of several subsystems (see [18-21]). First a mass 
flow-metered gas mixing system delivers a fuel mixture (sometimes N 2 -diluted), containing one of the 
gaseous hydrocarbon fuels, to the upper speaker-driver and then down to the nozzle exit. The fuel 
mixture enters a shallow cone-shaped plenum (~ 6-cm high at the center), bounded at the top by a 20- 
cm diameter polypropylene-coated speaker-driver diaphragm, and below by a machined metal plate. 
The plate has a 2.5-cm diameter hole, centered on the speaker axis, which is O-ring sealed to a 
vertically-oriented 2.2-cm i.d. Pyrex tube that effectively extends 38.3-cm to the nozzle exit. Similarly, 
mass flow-metered air is delivered through an identical speaker-driver system at the bottom. Thus in 
each experiment, nearly-uniform laminar jets of fuel mixture and air are impinged through matched 
pairs of 7.2 mm convergent Pyrex nozzles, to form a free-floating flame (after ignition); this flame moves 
vertically with any flow imbalance (hence free floating), and is always kept centered. The nozzle gap is 
fixed at two exit diameters. The nozzle area contraction ratio is ~ 9:1 with convergence over ~1 tube 
diameter, and each nozzle is slightly "recurved" near the exit. The resultant plug-flow velocity field, 
which exhibits a shallow dish-shaped central depression (consistent with numerical predictions by 
Rolon et al.), has been previously characterized using LDV and PIV measurements [18,19]. The Pyrex 
nozzles, insulated by blocks of silica foam insulation, are mounted in a rigid ceramic fiber box with three 
Pyrex windows and a porous top of sintered metal or ceramic fiber. Argon or N 2 bath gas, dispersed 
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radially via a jet near the bottom of the box, prevents/inhibits extraneous combustion outside the central 
impingement region, and thus minimizes adverse flame attachment, buoyancy, and visibility effects. 

Fuel and air component flows were hand-controlled using micrometer valves, and monitored by the 
calibrated mass flowmeters. 

Steady State Flame Strength Measurements 

To obtain steady-state extinction of a CFDF, the mass flows of fuel mixture and air are slowly increased 
simultaneously, so the disk flame, located primarily on the airside, is always centered and free-floating 
[16-21]. Manual control of flame centering is monitored visually, and also through a TV display from a 
horizontally oriented focusing schlieren system [18]. Whenever N 2 -diluted fuels are used, the fuel is 
fixed at a target rate and N 2 is increased. In this case, because metered fuel and diluent flows are 
blended in a small glass bead mixer, and then pass through a substantial “dead volume” in the fuel 
speaker plenum and tube/nozzle system, diluent flow is increased very slowly so the drift in mixture 
composition has sufficient time to reach the nozzle exit. After each extinction (blowoff), or rupture of the 
(very) flat disk-shaped-flame structure (monitored via focusing schlieren), a residual ring-shaped flame 
quickly establishes at a stable location centered near the stagnation point, where the degree of mixing 
is high. Mass flow rates of each component are then recorded. 

The global nozzle-area-average jet velocity used to evaluate flame strength (FS) at extinction, l/ air , is 

calculated from the mass flow of dry “service” air standardized at 273 °K and 1 atm, and measured 
nozzle exit diameter, 7.2 mm. The L/ fue i is evaluated similarly from component mass flows. 
Corresponding Reynolds numbers are generally less than 1500, but are not so low that CFDF's are 
unacceptably thick or non one-dimensional. 

To restore the ring-like “tribrachial” or “edge” flame to a disk flame, respective flows are gradually 
decreased, so the slowly shrinking ring approaches ~1 jet diameter near the stagnation surface. At 
flame restoration the ring suddenly propagates inward and shifts axially (say, ~ 1 mm) to the airside. 
After flame restoration, another set of extinction/restoration measurements is obtained for replication 
purposes. 

In previous studies of H 2 — air and HC— air systems, flame restoration was found independent of jet 
diameter, and a large hysteresis in exit velocity existed between extinction and restoration [16-21]. It 
was concluded that flame restoration occurs as a velocity-limited piloted-reignition along a thin 
stagnated region containing inter-diffused jet flows; or, when the inward "stretched laminar burning 
velocity" finally exceeds the maximum outward radial velocity [21]. Recent very detailed numerical 
simulations [23,24,28] fully support our earlier simplified description. Note, for extinction-restoration 
hysteresis, the radial strain rate in the central stagnation region must always be smaller than that 
required for extinction. Otherwise blowoff and restoration will occur at essentially the same flow rates - 
which was observed for methane-air using (unnecessarily) small 2.7-mm tube-OJB’s [16]. 

Dynamic Flame Strength Measurements 

Dynamic extinction measurements were obtained using the same basic flow technique as in the steady- 
state measurements, except the twin speakers were driven in-phase (almost always, unless noted) by 
one of several pre-selected applied voltages, and frequencies. Thus, a waveform generator was first 
set to a desired frequency and a sinusoidal output of 1 volt amplitude. This became input to a variable 
speaker-amplifier. Twin analog outputs from the amplifier powered the twin speakers. The amplifier 
outputs were manually adjusted to selected peak-to-peak (pk/pk) voltages while viewing waveform 
displays on a LeCroy digital oscilloscope. Once the speaker-drivers were set identically, a dynamic 
Flame Strength measurement was obtained using the same flow procedures as discussed above. If 
flow reversal in the nozzles appeared likely, based on previous hot-wire flow calibrations (described 
below), and/or the flame was too unstable, based on visual and focusing schlieren observations and 
also difficulties in maintaining a centered and unanchored flame, lower pk/pk applied voltages were 
used. 
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During very early stages of the study (with B. Reid in 1998 [26]), dynamic extinction measurements 
were obtained for the 18% H 2 -air and CH 4 -air systems over a range of pk/pk speaker-driver voltages, 
at respective frequencies of 30, 60, 120 Hz; and later at 25, 50, 100 Hz. Thus some limited sets of 
dynamic flame strength data corresponded to known speaker-voltage waveforms. To analyze these FS 
data in terms of physically realistic air and fuel input velocities / strain-rates that effectively characterize 
extinction limits, independent velocity information at each nozzle exit was needed as a function of 
applied speaker voltage. (Note that limited sets of definitive measurements of dynamic axial strain rate 
near the airside flame edge have been reported by others [1 3-1 5].) 

Thus limited Hot-Wire velocimetry measurements of cold flows at the nozzle exits were obtained at 
three frequencies each, over two distinct time periods (2000 and 2001 ) [26]. In each case, pk/pk 
sinusoidal air and fuel cold-flow input velocity magnitudes were measured at the respective nozzle exits 
corresponding to various sets of pk/pk speaker voltage magnitudes, steady input mass flows, and 
applied frequency. The respective Hot-Wire velocity data sets were analyzed using three levels of 
empirical curve fits. Thus dynamic extinction data could be transformed and analyzed as functions of 
pk/pk sinusoidal air and fuel input velocity magnitudes. The velocity magnitude inputs were calculated 
as empirical functions of pk/pk speaker voltage magnitudes, mass flow rates of air and fuel, and 
frequency up to 120 Hz, using an expression for L/ a i r that produced nearly identical results compared 
with the later set of independently derived expressions for L/ a i r and Uf ue i (unpublished data, Clare Me 
Namara (2000), Princeton University; and Rachel Johnson (2001), The Air Force Academy). 

Upon combining the dynamic input velocity and extinction data for C 2 H 4 /N 2 -air systems, it was found, 
for any given raw (uncorrected) C 2 H 4 mass flow, i.e., 3, 5, or 8 SLPM, the measured flame strength, 

U a \ r , was a linear decreasing function of the pk/pk velocity magnitude of the sinusoidal input, pk/pk U a \ r . 
The resultant negative slopes of these data were primarily a function of frequency, but they also varied 
with the respective input mass flow rates of C 2 H 4 . (Later on, slope variations for the three different input 
mass flows of C 2 H 4 were effectively “normalized-out” using corresponding steady-state flame strengths 
of C 2 H 4 /N 2 -air, to produce a singular Dynamic Flame Weakening response for the C 2 H 4 /N 2 -air 
systems; see Results section.) 

Errors in extinction limits stemmed from various sources; and some are compensated-for internally. For 
example, according to the ideal gas law, absolute jet velocities at constant mass flow vary linearly with 
input temperature. Hence mass inflow measurements are unaffected by local variations in gas velocity 
caused by jet temperature variations, and resultant “density-weighted” extinction limits have been found 
independent of such linear variations in gas velocity up to ~ 600 K [21]. Thus potential data scatter due 
to temperature variation is effectively nulled, which is especially helpful when heated fuel jets are 
required [17]. However, variations in atmospheric pressure do affect the data slightly. These have not 
been routinely accounted-for until very recently [29], where U a \ r calculated at 300K is now routinely 
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multiplied by the ratio ((29.92 in Hg)/(Barometric Pressure, in Hg)) . Although the same nozzle-OJB 
was used throughout this study, on an absolute basis the calculated jet exit velocities are sensitive to 
(measured) jet diameter squared, and strain rates are sensitive to jet diameter cubed. This can be 
important when results are compared from different-size OJBs. 


Results 


Steady State Flame Strength 

Steady-state extinction results for C 2 H 4 /N 2 -air CFDF’s, shown in Fig. 2, exhibit typical data 
reproducibility and variation of FS with input mole fraction, for a N 2 -diluted gaseous hydrocarbon fuel. A 
minor extrapolation of the asymptotic data trend, using a simple polynomial fit (previously well- 
established [1 6-21 ]), indicates a FS of 31 0 cm/s for 1 00% C 2 H 4 -air flames (at 273 K). This 
corresponds to an Applied Stress Rate of 310/0.72 = 430 1/s, and a global axial strain rate at the airside 
edge of 2 x 430 = 860 1/s at 273 K, for an idealized 1-D CFDF with uniform axial velocity input [21]. 


5 



The asymptotic limit represents a chemical kinetic limited reaction rate for a CFDF that is unaffected by 
the diffusion rate and thermodynamic properties of the N 2 . The indicated raw C 2 H 4 inputs, in 
uncorrected Standard Liters per Minute (SLPM), do not contain a factor of 0.69 for ethylene mass flows 
(based on heat capacity) that is otherwise applied to all calculated quantities. 

The steady-state FS for 100% C 2 H 4 -air can be compared with previous measurements for 100% CH 4 - 
air and 18% H 2 -air CFDF's obtained using the same 7.2 mm OJB, and 100% H 2 -air CFDF's that were 
characterized earlier using a series of scaled OJB’s [21]. Based on the previous data [21,26,27] and 
also an earlier FS of 1 1 3.6 cm/s @ 300 K for Chl 4 -air, the C 2 H 4 -air flame is 3.0 times stronger than the 
CH 4 -air flame. And based on the extensive H 2 -air nozzle-OJB results in [21], the C 2 H 4 -air flame is 
12.0 times weaker than the 100% H 2 -air CFDF. 

Dynamic FS, and Data Transformation using Limited Sets of Hot-Wire Velocimetry of Cold Flows 

Fig. 3 exemplifies a transformation / analysis of the dynamic extinction data, using the previous Hot 
Wire velocimetry data to quantify the effect of sinusoidal pk/pk exit velocity magnitudes (from sinusoidal 
pk/pk voltage magnitudes to speakers) on the dynamic Flame Strength of C 2 H 4 /N 2 -air CFDF's, for 3 
SLPM (raw) C 2 H 4 input flows. Thus the abscissa is based on our earlier [26] empirical analysis of Hot 
Wire data on nozzle exit flows, up to 120 Hz, and then an over-extrapolation of the empirical pk/pk exit 
velocity magnitude function to higher frequencies (beyond 120 Hz). 

The Fig. 3 data apply to a selected range of frequencies (100-500) over which the weakening effects on 
flame strength appear to grow large, and then become negligible. The data at both 100 and 200 Hz 
show the dynamic response of flame strength to “projected” pk/pk L/ air is linear, even though the 
magnitudes of oscillation are large compared with the mean L/ air . Checks at other flows and 
frequencies indicated similar linear responses, even under rare conditions where local flow reversal 
(which we carefully avoid) is known to have occurred (from HW data traces). The large change in 
projected dynamic response (slope) in Fig. 3, that occurs between 450 and 500 Hz, appears to suggest 
some sort of transition, where Flame Strength may no longer be affected. 

Thus the linear slopes in Fig. 3, which apply at particular frequencies and are well defined in numerous 
similar plots, are referenced to the magnitude of our earlier empirical pk/pk exit velocity function 
(inappropriately extrapolated to frequencies higher than 120 Hz), and seem to represent the sensitivity 
of dynamic flame strength to projected pk/pk sinusoidal velocity inputs at various frequencies. These 
slopes were initially thought to provide a useful measure of the dynamic response of flame extinction to 
imposed velocity oscillations, as a function of frequency (summarized later in Fig. 5). 

Fig. 4 shows a wide range of dynamic extinction results, as an extension of the steady state results in 
Fig. 1 . It represents ~ 200 C 2 H 4 /N 2 -air CFDF measurements made early in the program. These were 
obtained using various sinusoidal speaker-voltage inputs that produced various predicted axial velocity 
input magnitudes, pk/pk L/ air , at preset frequencies. For certain mid-range frequencies, the data 
illustrate significant impacts of oscillatory velocity inputs on dynamic Flame Strength (up to a factor of 
2), over a wide range of raw C 2 H 4 inputs and resultant C 2 H 4 /N 2 mole fractions. At frequencies of 500 to 
1600 Hz there is very little effect of input oscillation on flame strength compared to 0 Hz (steady-state). 
At very low frequencies (e.g. 4 and 10 Hz) some minor “apparent” weakening is evident, but this 
occurred at relatively large amplitudes of oscillation that added to the imposed mean velocity (strain 
rate). 

Fig. 5 summarizes the resultant “projected dynamic velocity response” of Flame Strength, obtained 
from the 200 early measurements of C 2 H 4 /N 2 -air extinction in Fig. 4, and based on the empirically 
derived velocity magnitude projections from Hot Wire data. Because the respective data sets for 3, 5, 
and 8 SLPM C 2 H 4 exhibit very similar trends and slightly more scatter (shown later), data averaging is 
used to define the apparent effect of frequency more clearly. The 4 and 1 0 Hz, data are believed to be 
essentially quasi-steady, so that potentially limiting molecular / free-radical transport processes are fast 
enough to follow the oscillations in input strain rate. In such cases, flame extinction should occur near 
peak input velocity, which approximately equals U a \ r + (pk/pk L/ air )/2. At higher frequencies, however, 
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the “dynamic velocity response” of Flame Strength grows more negative until, just beyond 250 Hz, it 
appears to decrease steeply by a factor of four. And shortly after this fall-off, a very steep recovery 
occurs between 450 and 500 Hz, which appears to define a highly weakened “notch” region terminated 
by a steep high-frequency cut-off limit. The reproducibility of these extinction data is remarkably good 
at fixed frequencies (typically + 3%), even at highly sensitive frequencies, e.g. 475 Hz. At frequencies 
above 500 Hz, the CFDF's become minimally responsive, and tend asymptotically to become totally 
unresponsive with respect to extinction. 

At this point, the likelihood of significant acoustic resonance effects became clear. The small waviness 
in the spline fits of Fig. 5 suggested minor resonance phenomena occurred over some regions; 
however, the large and abrupt changes in dynamic velocity response around the “notch” region 
suggested that major resonance phenomena were also occurring. A simple calculation pointed to the 
problem. The 38-cm length of Pyrex tube to each nozzle exit, plus ~ 6-cm of dead space between the 
tube and the speaker diaphragm in the plenum region, coincides with two wavelengths of sound at 395 
Hz (= (348 m/s)/(0.44 m)(1/2), which would promote constructive interference within each speaker-tube 
system. Although obvious anomalous changes in sound intensity were not recalled when frequency 
was changed at fixed applied voltage, it was clearly important to conduct a sound measurement survey. 

Probe Microphone Measurements 

A wide-frequency-range calibrated Probe Microphone system was used to measure the localized sound 
pressure field, at the air (and fuel) nozzle exits, as a function of frequency, both without flows, and 
later with fixed flows of 5 SLPM N 2 . The purpose of the initial PM measurements at zero flow [26] 
was to characterize resonance phenomena and provide a quantitative basis for characterizing dynamic 
flame response. Later, it was recognized that the dynamics of each speaker diaphragm, and resultant 
sound pressure measurements, were significantly sensitive to the pressures of gas flowing through 
each speaker plenum to the respective nozzles located downstream. Notably, this was already known 
to be true for the earlier Hot Wire velocity calibration measurements. 

Both sets of the Fig. 6 results show substantial effects of acoustic resonance in the system, equivalent 
to constructive interference for 2, 1 , 1/2, and 1/3 wavelengths. And, Fig. 6 also shows that significant 
acoustic pressure differences exist, at both low and middle frequencies, between PM data obtained with 
zero flows and with fixed flows of 5 SLPM N 2 . More detailed tests (not shown) at key frequencies, and 
with flows of N 2 between 0 and 8 SLPM, indicated that most of the acoustic pressure differences 
occurred from zero to relatively low flows, and the detailed data obtained at 5 SLPM adequately 
represented the experimental range of flows encountered. 

Finally, while it is recognized that the Euler equation can be solved numerically with appropriate 
boundary conditions, to obtain the corresponding oscillating exit velocity field in our confined system, 
where acoustic impedance is important, this has not been attempted. Meanwhile, the dynamic 
extinction results are analyzed in a similar manner as before, but now using the 5 SLPM sound 
pressure magnitude data as follows, instead of the frequency-limited velocity magnitudes. 

Transformation of Dynamic Data using Probe Microphone Measurements 

Fig. 7 shows the same dynamic extinction data as in Fig. 3, except the new abscissa is measured pk/pk 
sinusoidal pressure magnitude at the air nozzle exit, based on zero cold flows. The data are similarly 
linear at each frequency, and thus slopes denoted as “dynamic pressure response” are defined. Since 
no evidence of significant non-linearity was found in the complete dynamic extinction data set for 
ethylene, all the data can be analyzed (as before) as a function of frequency. Although corresponding 
pk/pk sinusoidal velocity magnitudes are clearly preferable to pk/pk pressure magnitudes, to allow a 
more physically meaningful analysis of the results, this would require a solution of the Euler equation for 
the OOJB system, which is presently unavailable. 

As a final step, the same dynamic extinction data shown in Figs. 3, and then reanalyzed in Fig. 7, are 
corrected in Fig. 8 as a function of pk/pk sinusoidal pressure magnitude at the air nozzle exit, based on 
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5 SLPM flows. The resultant ’’corrected” dynamic pressure response slopes in Fig. 8 (for 3 SLPM 
C 2 H 4 ) were then coupled with numerous other slope results from similar data plots (not shown), and 
finally assessed as a function of frequency. As with the first analysis using PM data based on zero cold 
flows, the original “notch” in the dynamic-velocity-response data (Fig. 5) continued to be entirely absent 
as a result of normalizing the dynamic extinction data using PM data based on 5 SLPM flows. At this 
point the entire set of dynamic pressure response results for the C 2 H 4 /N 2 -air system, based on 3, 5, 
and 8 SLPM C 2 H 4 flows, could be examined in an attempt to rationalize the results. 

Refinement of the Flame Response Function for the C 2 H 4 /N 2 -Air System 

Once the entire set of dynamic Flame Strength data for the C 2 H 4 /N 2 -air system was reanalyzed, to 
obtain linear slopes as a function of pk/pk pressure magnitude based on the 5 SLPM PM calibration, a 
plot of the raw Dynamic Flame Weakening data could be examined as a function of frequency, as 
shown in Fig. 9. Several important features are evident. First the respective raw DFW data for the 
three different ethylene flow rates exhibit significant differences, especially in the 8 to 50 Hz range. 

This suggests the results might be normalized by Steady-State Flame Strength to collapse the data to a 
singular function (discussed below). Second, a fairly well defined region of Quasi-Steady responses 
appears to exist at frequencies up to 20 Hz. Third, a much different mid-frequency transition is evident, 
and the “high-frequency” cutoff beyond 200 Hz is followed by an asymptotic approach to “complete 
insensitivity” at frequencies from ~ 300 Hz to 1600 Hz, the limit of this study. And, as noted above, the 
former “notch” in dynamic-velocity-response data effectively disappeared entirely as a result of 
normalizing the dynamic Flame Strength data with PM pressure data. 

Fig. 10 shows a complete set of fully normalized DFW results for the ethylene-air system. These were 
obtained after the raw data in Fig. 9 were normalized by Steady-State Flame Strengths for the 
respective ethylene flows of 3, 5, and 8 SLPM (U air ss = 210.8, 256.4, and 289.6 cm/s); the results were 
also multiplied by 1 00 to define the percent Dynamic Flame Weakening per unit of pk/pk pressure 
magnitude (Pa). Thus the raw Dynamic Flame Weakening results in Fig. 9 collapsed to a single curve 
in Fig. 10, which characterizes the ethylene response independent of its mole fraction in mixtures with 
N 2 . Note the DFW data in the 8 to 100 Hz range were most affected, with a substantial reduction of 
“data scatter” for the 3, 5, and 8 SLPM ethylene mass flows. These data represent a respectably large 
range of Steady State Flame Strengths and ethylene mole fractions. 

Normalized Dynamic Flame Weakening of Methane-, Ethane-, and Propane-Air CFDF’s 

Fig. 1 1 shows a collective plot of the entire set of normalized Dynamic Flame Weakening results for 
hydrocarbons in this study, including recent data for the three additional systems. Clearly, the collective 
results show that methane flames are weakened to a much greater extent than ethylene flames. The 
propane and ethane flames exhibit intermediate behavior. Also, the “knee like” transition from quasi- 
steady extinction, to “progressive strengthening” of flames in the “falloff’ region, occurs at progressively 
higher frequencies as the fuel changes from methane to ethylene. Notably this transition was predicted 
by Egolfofopoulos and Campbell [11], who used an analogous approach to Stokes’ Second Problem 
that describes the attenuation of oscillations by diffusion. Because all the data are normalized with 
respect to Steady-State Flame Strengths (e.g., L/ airss = 1 14.3 cm/s at 300 K for methane), the absolute 
differences in Flame Response between ethylene and methane are somewhat smaller. 

Finally, Fig. 12 shows a summary crossplot of all the fully normalized Dynamic Flame Weakening vs 
frequency results in Fig. 11, as a function of Steady State Flame Strength, L/ air>ss , for all the HC-Air 
systems studied. Clearly, the respective DFW’s of CH 4 -air, C 3 H 8 -air, C 2 H 6 -air, and C 2 H 4 /N 2 -air flames 
(in order of decreasing DFW) correlate remarkably well with the Steady State FS’s (in order of 
increasing SS FS's). The resultant power law exponents vary from -1 .5 for Quasi-SS limits at < 1 5 Hz, 
to - 1 .0 at 100-200 Hz, to ~ 0 beyond 300 Hz. Note that the Quasi-SS limits place an upper bound on 
Dynamic Flame Weakening. Although close inspection of the propane-air results in Fig. 11 indicates 
the relative decrease in DFW is clearly greater in the / range, 15-80 Hz, than for the other fuels, the Fig. 
12 correlation still shows that Dynamic Flame Weakening of all the flames in this study, over the entire 
range of important frequencies, is always greatest for flames that are weakest in the Steady-State. 
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Finally, the above departure of propane results may be a first indicator of additional chemical kinetic 
effects, and/or it may reflect a change in CFDF structure in which the flame zone is closer to the 
stagnation plane. 

Effects of Input Phase on Dynamic Flame Response 

It was considered important, from both theoretical and practical standpoints, to measure the possible 
effects of fuel-air input phase differences on dynamic Flame Response. Note that an earlier single 
waveform preamplifier with twin output amplifiers was replaced by a system of twin preamplifiers, and 
an amplifier with a settable phase lag. In all cases the respective voltage waveforms applied to the twin 
speakers were monitored using a digital oscilloscope (LeCroy) system. 

Fig. 13, which shows 3 SLPM ethylene results obtained at 150 Hz, indicates there is a mild phase angle 
effect over the very large range investigated. Each data point represents at least duplicate 
determinations, and at zero phase the data point represents a grand average of measurements 
throughout the study. 

Fig. 14 shows the additional phase-effect data obtained, which includes frequencies of 50, 75, 100, 150 
and 200 Hz. Although the data trends for the two higher frequencies are nearly flat (including the 1 50 
Hz data discussed above), the average trends for the lower three frequencies trend upward and show 
increasing scatter, especially the 75 and 50 Hz results. The apparent upward trends between zero 
phase and 180 degrees might be rationalized on the basis that, e.g., a 180 degree shift at 50 Hz 
represents a “push-pull” oscillation that is effectively 100 Hz. Despite the data scatter, best-fit zero- 
phase intercepts agree very well with corresponding averaged values in the previously evaluated zero- 
phase data for 3 SLPM ethylene, shown in Fig. 10 and 11. 

Implications for Scramjet Combustors 

Possible implications of the dynamic results for scramjet and other high speed combustion systems are 
certainly varied and difficult to determine, considering the known complexity of flameholding processes, 
reviewed in [22]. One relatively obvious approach for scramjet flameholding is to consider fundamental 
resonant frequencies in cavity pressure oscillations caused by ducted supersonic flows, and resultant 
turbulent shear layer flows, over both open and closed cavities [25]. For example, application of a 
“modified Rossiter” expression, reviewed in [25], predicts a fundamental frequency of 2960 Hz for a 
typical free stream air velocity of 1500 m/s over an open cavity of 10-cm. length. Because this 
frequency is considerably beyond the present high-responsivity range for hydrocarbon-air combustion- 
extinction, such induced oscillations may not affect open-cavity-based flameholding in this idealized 
case. Nevertheless, careful acoustic measurements and more detailed analyses of flameholding in 
various other fueled cavities and subsonic wake flows, especially in the vicinity of possibly-oscillating 
fuel jets, are needed to fully assess the potential for flame weakening effects in both test engines and 
flight vehicles. 


Findings and Concluding Remarks 

A unique idealized experimental study of the subject fuel-air systems was performed using an 
Oscillatory-input Opposed Jet Burner (OOJB) system. The study adds two new hydrocarbon data sets, 
and includes a significantly refined data analysis, compared to the author’s previous paper, AIAA-2003- 
4634. Extensive dynamic-extinction “Flame Strength” (FS) measurements were obtained on 
unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's), stabilized by steady input 
flows and perturbed by superimposed in-phase sinusoidal velocity inputs at both the fuel and air nozzle 
exits. Both limited Hot-Wire (HW) measurements of velocity fluctuation magnitudes at the nozzle exits, 
and much more extensive and accurate Probe Microphone (PM) measurements of acoustic pressure 
magnitudes, were obtained with cold flows and then used to analyze the dynamic Flame Strength data. 
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It was very important to use PM measurements obtained when significant gas flows were maintained (a 
nominal mean flow of 5 SLPM was used), rather than earlier PM measurements obtained under zero- 
flow conditions, as was done in the author’s previous paper, AIAA-2003-4634. This is because each 
speaker diaphragm needs to be positively pressurized by flowing fuel and air, to replicate normal run 
conditions with speaker-imposed oscillations during extinction measurements. 

The steady state and dynamic CFDF extinction measurements in this study, coupled with the respective 
dynamic cold-flow nozzle-input calibrations, led to the following findings and conclusions: 

(1) The HC— air flame weakening results are unique but fundamentally consistent with limited 
published results on the structure and extinction (rare) of forced unsteady CFDF's. Available numerical 
simulations and sparse data indicate diffusive responses of key species (e.g. H-atom, H 2 , O and OH) 
control phase lags (relative to inputs) in concentration and temperature waves, which may weaken a 
flame. At high f, rate-limited diffusive responses ultimately lead to flame insensitivity, with no 
consequent weakening. 

( 2 ) The 0 Hz (Steady-State) Flame Strength for 100% C 2 H 4 -air can be compared with FS's for 100% 
CH 4 -air and 18% H 2 -air CFDF's obtained using the same 7.2 mm OJB, and 100% H 2 -air CFDF's 
characterized earlier using scaled OJB’s [Pellett et al, Comb. & Flame, 112:575-592 (1998)]. Based on 
the present FS of 1 14.3 cm/s at 300 K for CH 4 -air (2 U a JD n = 318 1/s, global axial strain rate), the 
C 2 H 4 -air flame is 2.98 times stronger. And yet, based on previous (and more recent) H 2 -air results, the 
C 2 H 4 -air flame is 12.0 times weaker than a 100% H 2 -airCFDF. 

( 3 ) Dynamic FS varied linearly with pk/pk U air (velocity magnitude, HW), and also pk/pk P (pressure 
magnitude, PM) at nozzle exits. A complete set of PM cold-flow measurements effectively normalized- 
out (accounted-for) large acoustic resonance apparatus effects in the dynamic FS data. The Dynamic 
FS data were also normalized by 0 Hz (SS) FS. Thus Dynamic Flame Weakening (DFW) was defined 
as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = -100 d(L/ air / U a „ i0Hz ) / d(pkpk P). 

( 4 ) For the C 2 H 4 /N 2 -air system: From 8 to 20 Hz, DFW is constant at 8.52 + 0.2 (% weakening)/Pa, 
and surprisingly independent of C 2 H 4 /N 2 mole fraction. This reflects a “quasi-steady flame response” to 
an acoustically induced dU a „ldP. From 20 to ~ 150 Hz, the C 2 H 4 /N 2 -air flames weaken progressively 
less, with an inflection at ~ 70 Hz, and become asymptotically insensitive (DFW ~ 0) at ~ 300 Hz, which 
continues out to 1600 Hz. The DFW of CH 4 -air flames follows a similar pattern, but shows much 
greater weakening than C 2 H 4 /N 2 -air flames; i.e., the quasi-steady DFW of CH 4 -air (8 to 15 Hz) is 44.3 
%/Pa, or ~ 5x larger, even though the 0 Hz steady-state FS is only 3.0 x smaller. The quasi-steady 
DFW's of C 3 H 8 -air and C 2 H 6 -air are intermediate at 34.8 and 20.9 %/Pa, respectively. 

( 5 ) Thus DFW responses for f< 15 Hz are effectively quasi-steady. At higher frequencies, DFW 
remains significant but decreases progressively, apparently reflecting increasing internal phase lags in 
diffusive transport. A practical high-f limit for DFW is effectively ~ 200 Hz, beyond which the flames 
become essentially insensitive to input flow oscillations out to 1600 Hz, the limit of the study. 

(6) Respective DFW's at various frequencies, of CH 4 -air, C 3 H 8 -air, C 2 H 6 -air, and C 2 H 4 /N 2 -air flames 
(in order of decreasing DFW), correlate remarkably well with Steady State FS's (increasing SS FS's). 
Power law exponents vary from -1 .5 for Quasi-SS limits at < 1 5 Hz, to - 1 .0 at 1 00-200 Hz, to ~ 0 
beyond 300 Hz. The Quasi-SS boundary marks the greatest DFW that can be observed. 

( 7 ) The DFW of C 2 H 4 /N 2 -air is independent of N 2 dilution over a relatively wide range, due to 
normalization of the (PM-normalized) dynamic Flame Strength by SS Flame Strength. This is likely true 
for the other hydrocarbon-air systems. Thus the characteristic Dynamic Flame Weakening of each 
hydrocarbon-air system appears unique to the combustion kinetics of that system. 

(8) Finally, systematic Fuel-Air input phase effects were measured for the C 2 H 4 /N 2 -air system, but the 
relative changes in DFW appeared both complex and relatively small. 
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In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and 
correlate closely, but non-linearly, with Steady-State Flame Strengths. The Dynamic Weakening of 
methane flames is much greater than for ethylene flames, up to nearly 200 Flz. For reasons 
unknown, the Dynamic Weakening of flames didn’t increase at intermediate frequencies (e.g., at 20- 
50 Flz) compared to low frequencies (< 1 5 Flz), where quasi-steady weakening appears to dominate. 
Quasi-steady weakening ostensibly corresponds to the peak input strain rate just exceeding the steady- 
state extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in 
the fall-off region. 

The fundamental results of this study may apply to airbreathing scramjet engines in which critical early 
stages of localized subsonic flameholding are associated with fuel injection processes that “feed” 
recirculation zones having “sufficient” residence-time distribution. Important phenomena may include 
acoustic resonance in cavity flameholders and fuel-injector wakes; near-field coupling with subsonic 
vortex shedding; and dynamic coupling with unsteady fuel injection systems. Incipient flames may be 
weakened or extinguished by acoustically driven strain-rate and frequency sensitive processes, such as 
Fl-atom diffusion in a “laminar-flamelet-like” reaction zone. Because such limitations may cause the 
loss of “robust” flameholding, and possibly generalized flameout, it appears important to characterize 
these dynamic effects and assess the potential for occurrence in scramjet and other combustor 
designs. Measurements of acoustic fields are needed to assess possible effects of localized acoustic 
fields on critical flameholding in ground-based and flight tests of scramjet engine configurations. 
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Fig. 1 . Schematic of Oscillatory Opposed Jet Burner (OOJB) system with 
twin 20-cm speaker-drivers. Diode laser system is passive in this study. 
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Fig. 2. Flame Strength at Steady-State extinction of C 2 H 4 /N 2 -Air Counterflow 
Diffusion Flames (CFDF's), using 7.2 mm Pyrex Nozzle OOJB. U a ir is 
standardized at 273 K, 1 atm. 
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Fig. 3. Dynamic extinction of C 2 H 4 /N 2 -Air CFDF’s using 7.2 mm Pyrex nozzle OOJB, with axially-applied in-phase sinusoidal 
velocity inputs. Abscissa is based on an empirical analysis of Hot-Wire Velocimetry data on nozzle-exit cold-air flows, from 
25-120 Hz, and then an extrapolation of the tested empirical pk/pk exit velocity expression to higher frequencies. 
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Fig. 4. Steady-State and dynamic extinction of C 2 H 4 /N 2 -Air CFDFs, using 7.2 mm Nozzle 
OOJB with 1-10 SLPM C 2 H 4 Flows, and axially-applied in-phase sinusoidal velocity inputs of 
C 2 H 4 /N 2 and air at varied amplitude and frequency. Early data up to 1 1/02 are included here. 
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Fig. 5. Dynamic Velocity Response for extinction of C 2 H 4 /N 2 -Air CFDF’s, using 7.2 mm Pyrex nozzle OOJB 
with 3, 5, 8 SLPM C 2 H 4 flows, and axial sinusoidal velocity inputs. Dynamic Velocity Response is based on 
pk/pk nozzle exit air velocities that were originally derived from Hot Wire measurements on cold air and N 2 
flows at 25-120 Hz, and then empirically projected at higher frequencies. 
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Fig. 6. Normalized Probe Microphone response of sinusoidal speaker-driver sound pressure at exit plane 
(r/a = 0.5) of 7.2 mm Pyrex (air) nozzle, for both 0 and 5 SLPM Air and Fuel flows (air also used for fuel 
side), using calibrated Probe Microphone with both 50 and 100 mm long probe tubes. The pk/pk voltage 
magnitude applied to the speaker-drivers is normalized to 1 .0 volt in this calibration. 
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pk/pk Sinusoidal Pressure at Nozzle Exit, Pa 


Fig. 7. Dynamic Extinction of C 2 H 4 /N 2 -Air CFDF’s, for 3 SLPM C 2 FI 4 (raw mass flows), using 7.2 
mm Nozzle OOJB with axially-applied sinusoidal velocity inputs. Abscissa is based on Probe 
Microphone pressure measurements of sound at nozzle exit with zero mean flow, shown in Fig. 6. 
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Fig. 8. Dynamic Extinction of C2F4/N2-Air CFDF’s, for 3 SLPM C 2 H 4 inputs, using 7.2 mm 
Nozzle OOJB with axially-applied sinusoidal velocity inputs. Abscissa is based on Probe 
Microphone pressure measurements of sound at nozzle exit with 5 SLPM flow, shown in Fig. 6. 
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Fig. 9 . Raw Dynamic Flame Weakening for extinction of C2H4/N2 vs Air CFDFs, using 
7.2 mm Pyrex Nozzle OOJB, referenced to axially-applied sinusoidal velocity (pressure) 
Inputs probed by Microphone (cal. @ 5 SLPM) 
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Fig. 10. Fully normalized Dynamic Flame Weakening of C 2 H 4 /N 2 -Air CFDF’s in Fig. 9, using 7.2 
mm Pyrex Nozzle OOJB and measurements of axially-applied sinusoidal velocity inputs probed by 
Microphone with 5 SLPM flows. Note that normalization collapsed the three sets of ethylene results 
to an equivalent response that is independent of dilution by N 2 . 
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Fig. 1 1 . Dynamic Flame Weakening of all HC— Air CFDF systems studied: C 2 H 4 /N 2 , C 2 FI 6 , 
C 3 H 8 , and CFI4 vs Air, using 7.2 mm Pyrex Nozzle OOJB and measurements of axially- 
applied sinusoidal velocity inputs probed by Microphone with 5 SLPM flows. Note the three 
sets of ethylene results are equivalent, and independent of dilution by N 2 . 
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Fig. 12. Summary crossplot of Dynamic Flame Weakening vs frequency results in Fig. 11, as a 
function of Steady State Flame Strength, for all FIC-Air CFDF systems studied. All results are 
referenced to axially applied sinusoidal velocity inputs probed by microphone with 5 SLPM 
flows. The Quasi-Steady-State results for f < 15 Hz represent an upper bound on DFW. 
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Fig. 13. Effect of Fuel-Air Phase difference on dynamic extinction 
response of C2F4-Air CFDFs, using 7.2 mm Pyrex Nozzle OOJB with 3 
SLPM C 2 FI 4 Flows, and 150 Flz axially applied sinusoidal velocity inputs 
probed by Microphone. 
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Fig. 14. Effect of Fuel-Air Phase difference on dynamic extinction Flame responses of 
C 2 Fl 4 -Air CFDFs at four key frequencies, using 7.2 mm Pyrex Nozzle OOJB with 3 
SLPM C2H4 flows, and axially-applied sinusoidal velocity Inputs probed by Microphone. 
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